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Abstract 



A theory of the Kondo effect in quantum dots at zero temperature in the 

presence of arbitrarily intense AC potentials is presented. We generalize the 

Friedel-Langreth sum rule to take care of charge conservation and propose a 

consistent procedure to study a time dependent Anderson Hamiltonian. The 

effect of the AC potential on both the quantum dot density of states and the 

linear conductance shows the importance of using a theory which describes 

intradot finite interaction and nonperturbative effects. 
PACS numbers: 73.40.Gk, 72.15Qm, 85.30Vw, 73.50.Mx 
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The competition between electron-electron interaction and quantum-mechanical hy- 
bridization between the delocalized electrons in a non-magnetic metal and the unpaired 
electrons of a magnetic impurity leads to the Kondo effect [[J. It has been predicted |2|-|9| 
that the transport at low temperatures through a single level coupled to two reservoirs 
containing Fermi liquids |1U|| , is governed by a Kondo-like singularity existing in the quasi- 
particle density of states (DOS) at the Fermi level €f of the reservoirs. This has been 
experimentally performed recently |TT| , [T2|| by means of a quantum dot (QD) coupled to two 



leads by tunneling barriers. When i) the QD is small to have well separated levels, ii) the 
tunneling rate is such that the broadening of QD levels is typically one order of magnitude 
smaller than the QD charging energy and in) the number of electrons within the QD is odd, 
then, there is an unpaired electron which is free to form a singlet state with the electrons 
in the leads. This is the physics behind of the Kondo effect that is well described by the 
low energy excitations of the Anderson Hamiltonian for the unpaired electron while all the 
other electrons are effectively described by a potential of the QD coupled to the leads. The 
excitations around ep are quasiparticles having a lifetime determined only by single particle 
properties without any many-body contribution. One electron at €p, becomes scattered by 
the dot suffering a phase shift which is proportional to the exact displaced charge 5n. As 
a consequence of this Friedel-Langreth (FL) theorem fl3| the conductance takes the value 
Q = sin 2 (n5n)2e 2 /h 0,f|. The interest of the Kondo physics in QD resides in the possibility 
of externally controlling parameters as the on-site interaction U, the QD level 6q or the 



broadening T of the DOS, due to the dot-leads coupling fTl]J12|| . It is also possible to study 
the Kondo singularity itself by applying a bias in order to mimic the QD DOS by the dl/dV 
curve 0,0,0. 



In this letter we address the following question: since the quasiparticles with large lifetime 
are those at ep, what happens to the Kondo problem at zero temperature when there is a mix- 
ing of states with energies differing from ep by a finite amount? In the case of QD, this can be 
actually made by means of an AC gate voltage [pl|-|r7|]. The problem is very interesting and 
far from trivial because the frequency range experimentally accessible (from ~ 100 MHz to 20 



GHz) |1| is roughly that of the Kondo temperature Tk = VUTexp{— n[\eF— £o\(U +e )]/TU} 
(from few mK to ~ IK) [0]. 

The aim of our work is to describe the QD DOS within a framework where charge 
conservation is taken into account in order to describe properly the Kondo singularity. 
Therefore, our first task is to generalize the FL sum rule to the case of a time dependent 
situation. The Anderson Hamiltonian is solvable through Bethe ansatz [19| or Quantum 



Monte Carlo methods pC§ but a reliable and simple method to obtain dynamical properties 
at low temperatures in the whole range of interactions (U/T) is not available. Previous efforts 
concentrated in the U — > oo limit where a Non-Crossing Approximation |jj |l4 , 17| can be 



made for high and intermediate temperatures. However, such approximations do not give the 
exact local Fermi-liquid properties when T — > and cannot describe the transition from the 
weak-correlation to the strong-correlation regime. On the other hand, finite U perturbation 
theory P, |2T| , |2~2"|| describes properly the symmetric case but presents clear anomalies away 
from this special situation. Our method is an extension of those |7|,p3| developed in the static 
case where an effective selfenergy eliminates consistently the problems of the perturbative 
approaches. 

The time dependent Anderson Hamiltonian reads: 

H = e kC ] Ku c k)(J + ( e * + V AC cosuj t) d\d„ 

k£{L,R},(T o" 

+ E V k {cl a d a + dlc k ^) + Ud\d^d\d l . (1) 

ke{L,R},a 

d\ creates an electron with spin a in the dot while c k a creates it in the lead with energy t k 
(k labels the rest of quantum numbers). The AC voltage with intensity Vac and frequency 
Uq modulates in time the relative position of the QD level e a with respect to ep. An 
eventual breakdown of the spin degeneracy would be represented by e CT ^ e_ CT . The coupling 
Vfc between the QD and the leads produces a broadening = —2Im[T,^ R \e + iS)] = 

27T Y^keL(R) | Vfc 1 2 (5(e — efc). The dc linear conductance at zero temperature deduced from a 



nonequilibrium Keldysh technique is [24 
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G= h T(e F ) (2) 

where T = T L + T R . Q is related to the time averaged DOS at the QD p CT (e) = —Im(A rT (e))/7r 
obtained from 

(M*)) = £-f*dt f dt.GKt^y^ (3) 

Z7T JO J-oo 

ti) being the QD retarded Green's function. Coherent tunneling of one electron at 
ep through the QD with AC produces a phase shift of its wave function which is related, 
according to the FL rule, to the displaced charge for each spin which is related with the QD 
occupation ((n a )) and with the number of lead electrons with ((n l a )) and without ((n l 0a )) 
coupling to the QD: 

Va (e F ) = 7r((n ff ) + (n l a ) - (n(, )(T ». (4) 
The phase shift is related with (^(e)) through the eigenvalues of the scattering T a matrix 
e 2^(e) = 1 + %{e) = 1 _ iJ^^^e)). (5) 



The displaced charge in Eq.(f|) can be obtained from the retarded Green's function as follows: 
i) = — ^ m {^k o-(^) 0) an d ( n o-) = ~I m (Ga(t>t))- ii) G K i s related to the retarded 
Green's function built up using G r a (t,ti) = G r a (t — txjexp —^Jl 1 dTVAc COSUJ o T > an d Hi) 
application of the operational rules given by Langreth |25| for Keldysh contour integration. 
This process, together with Eqs.(f|) and (|J), gives 

Imlni f; J 2 m ((3)G r a (e F - mKu ) \ = 

\m=—oo / 

-Im £ JttP) I" de T v( e) + !TT G "( e ~ mhuJ ^ 
_ J-oo lie — rrihujn) 



m=—oo 



( 00 IT/ I 2 \ 

where G r a (e) is the Fourier transform of the static retarded Green's function and J m the 
Bessel function of order m and argument (3 = Vac/^o- Eq.@ reduces to the usual FL rule 



T3fl for Vac ~ > 0. It establishes the charge neutrality constriction that must be verified by 
any approximation made to compute the retarded Green's function. This result implies a 
series of important consequences on the Kondo physics of the QD. 

Some approximations are made to compute the retarded Green's function 

^(6) = [6-6 (7 -S sp ( e )-S (T (6)]- 1 . (7) 

Correlation effects are included in the selfenergy S CT (e). In order to obtain S CT (e), we use 
an approximation which, in the static case [[7|,p3|, eliminates the pathologies of standard 
perturbation procedures. S CT (e) is calculated by an interpolation that gives correctly the 
limits U/Y — > and T/U — > and has good analytic properties both in e — > €f and 
e^ioo [0,0: 

ay ' <no,- CT )(l— (no,- CT )) v ; 

x (l _ j(l-(n | S (2) (e) A 1 (g) 

V (no,- CT )(l-(no,- CT ))C/ / c \ ) j \ J 

(tio.o-) is the mean field (MF) QD occupation and S^ 2 ^(e) is the second order selfenergy 
used in previous finite U calculations P,[7|, pB| , F2T| , F2^| . \i a is a shift of the QD level to an 



effective value determined selfconsistently together with (n a ) with the constriction of charge 
neutrality given by Eq.(||). Once these parameters are determined, spectral and transport 
properties are calculated from G r a (e). 

Hereafter we consider Y L = T R = V, e a = e_ CT = e and ep = 0. All our calculations 
are performed with parameters in the range of the experiments already available |TT|JT2 



in 



order to look for predictions on possible future experiments. First of all we analyze the effect 
of the AC field on the time averaged DOS that can be measured by means of the dl/dV 
curve of the QD as reported in the static case flIT] , [T2|j . Fig. [1] shows the DOS for a case 
with e 7^ —U/2. Results obtained from both Eq.@ for S CT (e) (continuous lines) and for the 
selfenergy up to second order (dashed lines) are shown. The selfconsistency for the selfenergy 
enhances the particle-hole asymmetry in both the static (Fig. |T]a) and high frequency (Fig. 
|l]d) cases. The lowest frequency (Fig. |I|b) corresponds to an adiabatic modulation of the 



level (hu)Q < r) which quenches the Kondo effect on the time averaged DOS. The two small 
peaks of the DOS correspond to the MF (lateral peaks) solutions separated by the extreme 
values of the harmonic potential. Fig. |I|c corresponds to the case (3 = 2 where the symmetric 
shape is recovered due to a transfer of spectral weight from the Kondo peak to the MF peaks 
and their satellites. Fig. [IJd shows the richest configuration with a strong asymmetry. This 
is better observed in Fig. ^| where the evolution of the DOS with Vac is analyzed. In the 
case of Vac = 0.4, one observes three types of peaks: the ones labeled as a correspond to the 
Kondo peak and its AC replicas. Those labeled as b correspond to the low-energy MF peak 
for Vac — and its AC replicas. These peaks evolve to higher energies and become narrower 
with increasing (3. Finally, the peaks labeled as c correspond to the high-energy MF peak for 
Vac — and its AC replicas. Once again, the evolution with (3 is similar to that of peaks b. 
Moreover, with increasing intensity, we observe a decrease of both the height and the width of 
the peak at the Fermi level. These effects are a consequence of the appearance of increasing 
new peaks and the necessary conservation of the total spectral weight. The reduction of the 
height affects the conductance as discussed below while that of the width implies a change 
in Tk- The behavior of the spectrum is not trivial {e.g. occupations change with (3), and 
clearly different to that obtained in a single particle picture of AC-assisted tunneling, where 
only simple replicas of the static DOS would appear. 

The enormous changes in the DOS when the occupation varies affect strongly to the 
conductance. Therefore, the effect of the AC potential on Q can not be described now by 
a single-particle model of AC assisted tunneling through a discrete state ]26|]. Fig [3| shows 
the AC effects on Q. At zero temperature and in the absence of AC potential [|§-§| Q has 
a broad peak of width given by U centered around eo = —U/2 and reaches the maximum 
value of 2e 2 /h PTfl . Fig ^| shows a reduction of the central peak (still centered around 
e = —U/2) and the appearance of symmetrically located satellites due to the AC. The 
position of the conductance satellites is inferred from the following argument: the DOS (left 
inset of Fig |3|) has a peak when the renormalized level of the QD is equal to Hujq. It has 
an AC satellite at (dashed line). This gives a large contribution to the conductance. For 
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low (3, the level renormalization is small so that the peak in Q is located at eo — hu>o — 0. 
Since the conductance is symmetric with respect to —U/2, that satellite has a companion 
at eo — — U — Tiujq. So, the separation between satellites is U + 2Tiujq as observed in Fig. 
|3] for Vac = 0.2. When the AC intensity increases, the QD level renormalization increases 
and the reasoning, still exactly valid for the renormalized level, only gives approximately the 
separation in eo between satellites. In the same range of values for U and T, the results are 
similar for Tiuoq = U/2, U/4 and 3 [7/20. It must be stressed that whereas the satellites in Q 
mark the AC-assisted transition from the empty orbital regime to the mixed valence regime, 
(inset at the left side of Fig |3|) the main broad peak corresponds to tunneling in the Kondo 
regime (inset at the right side of Fig|3|). This behavior in the conductance is different from 
the one predicted by a MF theory, where satellites at — Tiuj$ and hu — U also appear. The 
lack of these peaks is due to the different nature of the DOS when passing from the Kondo 
regime to the empty orbital regime (single peaked DOS centered at eo ) through the mixed 
valence one where the DOS has a single peak centered in a strongly renormalized position 
PIJTPj] . Our general result which depends on U contrasts with the results of Ref . |F| . They 



used a model with U — > oo and found peaks of half-width T in the conductance separated 
by hujQ, similar to previously reported MF AC-assisted tunneling. Another very interesting 
result drawn from Fig || is the non linear reduction of the central peak of the conductance 
when Vac increases. Since the whole behavior of Q with u is rather different to the linear 
dependence expected from a MF model, our result constitutes a proof of the importance of 
a finite U charge conserving theory for the AC Kondo effect. 

Finally, we pay attention to the effect of the AC potential on the occupation of the QD 
level. We show in Fig |^ the charge in the dot in terms of /3. Results are given for two different 
values of eo corresponding to initial occupations clearly above and below 1/2 (per spin) 
respectively. For increasing (3, (n a ) tends to 1/2 in a damped oscillatory way. For (3 such 
that Jo{(3) — > 0, the peak at eo in the DOS is quenched and its spectral weight transferred 
to the AC satellites. With the chosen parameters half of the satellites are above/below zero 
(with the same intensities), so that the QD occupation tends to 1/2. 



In conclusion, we have studied the AC Kondo effect in QD. We extend the FL sum rule 
to consider finite U as well as charge conservation. Our main findings are: i) changes in 
the shape and position of the peaks in the DOS in the presence of the AC potential, ii) 
variation of the DOS at ep and, consequently, of the linear conductance which has new 
features with respect to the usual MF theory. These are finite-U many-body effects. With 
regard to the possibility of experimentally observing the features described here, narrowing 
of peaks at the DOS due to AC effects are going to reduce Tr-. However, there is still 
margin for the measurements: for instance, in the mixed valence regime, T% ~ V ~ IK 
in available experiments Jl2] (without AC) while temperatures significantly lower can be 
routinely achieved. 

We acknowledge A. P. Jauho, L. Glazman, L.P. Kouwenhoven, A. Martin-Rodero and 
J.J. Palacios for helpful discussions. Work supported in part by the MEC of Spain under 
contracts PB96-0875 and PB96-0085 and by the EU via contract FMRX-CT98-0180. 
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FIGURES 

FIG. 1. Time averaged DOS with Vac = 0.2 for the asymmetric configuration U = 0.2, 
e = -0.15 and T = 0.025. a) without AC, b) hu = 0.01, c) huj = 0.1, d) huj = 0.35. 
Solid lines are calculated with the interpolation procedure for the selfenergy X CT (e) (Eq. ( ||)) and 
dotted lines with the selfenergy up to second order. 

FIG. 2. Time averaged DOS with hcoo = 0.35 for the asymmetric configuration U = 0.2, 
eo = —0.15 and T = 0.025 (as in Fig. [j]) and different intensities Vac- The labels a,b,c are 
discussed in the text. 

FIG. 3. Linear Conductance with Tiloq = 0.35, U = 0.2 and V = 0.025 for Vac = (continuous 
line), Vac = 0.2 (dotted line) and Vac — 0.4 (dashed line). The insets show the DOS with (dashed 
lines) and without (continuous lines) AC at the empty orbital (left side, eo = 0.35) and Kondo 
(right side, eo = —0.1) regimes. 

FIG. 4. Occupation (n a ) of the QD as a function of f3 for U = 0.2 and hu>o = 0.35. Results 
with the selfenergy up to second order (dashed lines) and with Eq.(||) (continuous lines) are shown 
for eo = —0.15 (curves above) and eo = —0.05 (curves below). 
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